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Abstract 

Schmutz Schaller's fV?\ , p. 201] conjecture regarding the lengths of the 
hexagonal versus the lengths of the square lattice is shown to be true. The 
proof uses results from (computational) prime number theory and from 
0- 

Using an identity due to Selberg, it is shown that the conjecture can 
in principle be resolved also without using computational prime number 
theory. By our approach, however, this would require a huge amount of 
computation, 

1 Introduction 

In [[H], p. 201] Schmutz Schaller, motivated by considerations from hyperbolic 
geometry, makes the conjecture that in dimensions 2 to 8 the best known lattice 
sphere packings have 'maximal lengths' and goes on to write: "In dimension 2 
the conjecture means in particular that the hexagonal lattice is 'better' than the 
square lattice. More precisely, let < hi < hi < ■ • • be the positive integers, 
listed in ascending order, which can be written as hi = x 2 + 3y 2 for integers x 
and y. Let < qi < qi < ■ ■ ■ be the positive integers, listed in ascending order, 
which can be written as qi = x 2 + y 2 for integers x and y. Then the conjecture 



is that qi < hi for i = 1, 2, 3, • • •." For an introduction to these matters see |TB 
For some progress regarding Schmutz Schaller's general conjecture in dimension 
2 see H (this case of the conjecture is also mentioned in [|], p. xxx]). 

For j > 1 let bj(n) = 1 if n is represented by the quadratic form X 2 +jY 2 and 
bj(n) = otherwise. The characteristic functions b\ and 63 are well understood. 
The following result was already known to Fermat. 

Lemma 1 A positive integer n is represented by the form X 2 + Y 2 if and only if 
every prime f actor p of n of the form p = 3 (mod 4) occurs to an even power. A 
positive integer n is represented by the form X 2 + 3Y 2 if and only if every prime 
factor p of n of the form p = 2 (mod 3) occurs to an even power. 
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(In general the natural numbers n that are represented by a quadratic form 
X 2 + mY 2 are rather more difficult to describe, cf. the beautiful book of D. Cox 
||.) Lemma [l] implies that hi and 63 are multiplicative functions. 

Let Bi(x) = J2 n < x bi{ n ) for i = 1 and i — 3. Schmutz Schaller's conjecture 
regarding the square versus the hexagonal lattice can be reformulated as follows 
in terms of B\ and B 3 . 

Conjecture 1 We have Bi(x) > B 3 (x) for every x. 

The first asymptotic result on Bi(x) goes back to Landau ||, who proved in 1908 
that 

B 1 (x)~C 1n -^=, (1) 



Vlogx' 



where 



C bl = ^= J] (l-P~ 2 r 1/2 = l II (1-P- 2 ) 1/2 - 0.764. (2) 

p=3(mod 4) p=l(mod 4) 

(Here and in the sequel the letter p is used to indicate primes.) Landau's proof 
uses contour integration. It is not difficult to use his method to show, cf. PD| , 
that for every k > 2 there exists constants C bl (2), ■ ■ • , C bl (k) such that 

B 1 (x) = C bl 1 + 1 h - — + O - — — . (3) 

Vlogx V. logx log fe 1 x \\og k xJJ 

This result can also be established by methods not using complex analysis, cf. 



T3, p. 288]. At the beginning of 1913 a then unknown Hindu clerk by the name 



of Ramanujan wrote in his first letter to Hardy @ that he could prove that 

B^) = C bl / -^= + 0(x^), (4) 

J 2 Vlogt 

for some e > 0. (For a reconstruction of Ramanujan's speculative argument see 
|l], PP- 60-66].) Note the similarity of Ramanujan's claim with the prime number 
theorem. From (|]) we infer that C bl (2) = 1/2 by partial integration. Shanks |2"T| 
showed, however, that C bl {2) ^ 1/2, thus disproving Ramanujan's claim. (Ra- 
manujan gave the correct formula and numerical approximation for C bl though.) 
The constants C bl and C bl (2) are known as the Landau-Ramanujan constant and 
the second order Landau-Ramanujan constant, respectively. For more on the 
evaluation of these constants see Section |^. For more on mathematical constants 
in general see e.g. 0. 

Ramanujan stated several claims similar to (|J) in his 'unpublished' manu- 
script on the partition and tau functions, see Section 6. All of them are disproved 
in | HJ| . It can be shown, however, that in each case Ramanujan's claims give the 
correct asymptotic main term. 

Similarly to (|3|) it can be shown that for arbitrary k > 2 there exist constants 
C& 3 (2), • • • , C b3 (k) such that 

ft(l) = cl . » ( 1 + ^P) + ... + ^w +0 ( i )), (5) 

Vlogx V log^ log x \log xj J 
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Table 1: -Bi(x) versus B 3 (x) 



where 

n (i- P - 2 r" 2 =^ n d-p- 2 ) i/2 « 0.639. 

V p=2(mod 3) p=l(mod 3) 

We thus arrive at the following conclusion. 
Proposition 1 Conjecture 1 is asymptotically true. 

Table 1 (copied from |22|] and verified by the second author) suggests that Con- 
jecture 1 is true for small x as well. The literature thus provides us with good 
indications that Conjecture 1 is true. The purpose of this paper is going beyond 
this and to prove that Conjecture 1 is indeed true. 

Theorem 1 We have B\(x) > B 3 (x) for every x, that is Schmutz Schaller's 
conjecture that the hexagonal lattice is 'better' than the square lattice is true. 



Landau's classical result (|l|) has been generalised in many directions, see [12] for a 



survey with over 50 references. Despite this rich history, nobody but the present 
author (in [|l(|) seems to have been concerned with proving effective results in 
this area, which is precisely what is needed to establish Theorem [l|. 



2 Preliminaries 

Let / be a multiplicative function from the natural numbers to lR>o- We define 

M f( X ) = Y.n<xf( n )i »f( X ) = T,n<xfi n )l n and X f( X ) = En<xf( n ) l °Sn. We 

denote the formal Dirichlet series Y^n=i f( n ) n ~ s associated to / by Lf(s). We 
define A/(n) by 



Lt(s) ^— ' n 

n=l 



n 



Notice that 

f(n)\ogn = J2f(d)Af(^). (6) 

d\n 

If / is the characteristic function of a subsemigroup of the natural integers with 
(1 <)gi < q2 < ■ ■ ■ as generators, then it can be shown that A/(n) = logg^ if n 
equals a positive power of a generator qi and Aj(n) = otherwise. For / = b\ we 
thus find, using Lemma [l], 

{21ogp if n = p 2r , r > 1 and p = 3(mod 4); 

logp if n = p r , r > 1 and p = l(mod 4) or p = 2; 

otherwise. 

For / = b 3 we find 

{21ogp if n = p 2r , r > 1 and p = 2(mod 3); 

logp if n = p r , r > 1 and p = l(mod 3) or p = 3; 

otherwise. 

From property (P) of A/(n), we easily infer that 

X f (x) = J2f( n )M^ (7) 

where ipf(x) = J2 n <x-^f( n )- ^ ne function Af and ^/ ar e analogues of, respec- 
tively, the von Mangoldt and the Chebyshev -^-function. 



3 Some related conjectures 

Unfortunately it seems that Mf is not a very 'natural' mathematical object, /if is 
on the other hand (as is amply demonstrated by browsing through the literature). 
For this reason we consider two additional conjectures: 

Conjecture 2 We have Xb 1 (x) > \b 3 (x) for x > 8. 



Conjecture 3 We have fJ>bi{ x ) > ^(x) for every x. 

Note that exp(A;, 1 (x)/2) is the product of all different lengths in the square lattice 
not exceeding yfx. Thus Conjecture 2 can be reformulated as stating that the 
product of the different distances not exceeding x occuring in the square lattice 
always exceeds the product of the different distances not exceeding x in the 
hexagonal lattice, provided that x > 2a/2. 

Conjecture 1 clearly implies Conjecture 3. Furthermore we have: 

Proposition 2 Conjecture 2 implies Conjecture 1. 

Proof. We have, for x > 2, 

M/W = r d 4M = m + rM^ it (8) 

/v J 2 _ log* logx J 2 tlog 2 t 
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Denote the latter integral by If(x). It is not difficult to show that h^x) > h 3 {x) 
for x < 8. Conjecture 2 then implies that the latter inequality holds for every x. 
The truth of Conjecture 2 together with @ then implies that Bi(x) > B 3 (x) for 
x > 8. By direct computation we then infer that the latter inequality holds for 
every x. □ 
Thus in order to establish Theorem [I], it suffices to establish Conjecture 2. 
From (0) and ipbi( x ) ~ x/2 as x tends to infinity it follows that Ab 4 (x) ~ /i^(a;)/2, 
as x tends to infinity An effective form of this relationship, together with an 
effective estimate for //^ (provided by Lemma H), then allows us to prove the 
main result of this paper: 

Theorem 2 The Conjectures 1, 2 and 3 are all true. 

Complications arise due to the fact that lim^oo B\ [x] / B 3 (x) = 1.1961377420 • • •, 
which is rather close to 1 and that ipbiiy) is not so close to y/2 for various ranges 
of small y (the convolutional nature of (0) forces us to take the small y range into 
account). 



4 The toolbox 

The following result from [[H]] will play a crucial role. It is in essence an effective 
version of Theorem A of 



Lemma 2 Let f be a multiplicative function from the natural numbers to M>o- 
Suppose that there exist constants D + and r, with r > 0, such that for every 

x > Xo, 

d-m*) < £ M ( E M5> _ T ,„ g I) < D+H{x) . (9) 

n<x I m<^ 

Then we have, for x > max{i , exp(D + )}, 



^log- ^ <^w<^i°g r ^ - ^ , (io) 

r 1 t 1 iix. 

log X log 2! 

where 

C f := \im (s-lYLAs). (11) 
In particular, if there exist constants C_ and C + such that 

A/(n) 



C_ < -rlogx < for x > 1, (12) 

n 

i/ien ( |T0| ) /ioWs true, /or x > exp(C + ), witt D_ = C_ and D + = C + . 
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Remark 1. From the proof of this lemma, fif(t)dt/t appears as a more easily 
estimated function than fif(x). Interestingly, Landau || in his proof of ([!]) using 
contour integration, estimates fj,t,i(t)dt/t rather than B\(x) itself. 
Remark 2. If \im x _ OQ (J2 n <x(-^f( n )/ n ~ T ^°S X )) exists, we denote this by Bf. 

Let us put 

t( n (log^-^+) T+1 ATT( n n] (logx-D_)- +1 

L(x,t,D_ } D + ) = — and U(x } t } D_ } D + ) = — . 

log x — D- log x — D + 

Thus we can write ( |10|) as CfL(x, r, D_, D + )/t < ///(-) < CfU(x, r, D_, D + )/r. 

Let r, s and c\ be given. At a few instances in the sequel we want to show that 
for every x > X2, with x 2 some explicit constant, we have fif(x/r) > ci/i g (x/s), 
where g satisfies the conditions of Lemma |2| with constants r, D'_ and D' + . By 
Lemma [| this leads us to consider inequalities of the form 

L(^, r, D>_, D+) > c 2 U(^, r, D'_, D' + ) (13) 

where all variables and constants are real numbers with r, r, s and c 2 positive, 
D- < D + , D'_ < D' + and x > xq := max{exp(_)+)s, exp(.D + )r}. We recall the 
following lemma from |TD| : 



Lemma 3 // logs + D'_ < D + + logr and Cp~3f) zs satisfied for some x\ > Xq, 
then (JT3) satisfied for every x > x\. If logs + D'_ > D + + logr ; and 

D'-D' \ D_-D + 

C2 1 + - ; + , : — < 1 



log(xi/s) — / log(xi/r) — D- 

for some x\ > xq, then (|T3|) is satisfied for every x > x%. 



We also need the following result about the difference between U(-, D + ) 
and L(f,±, D.,D + ). 

Lemma 4 Assume that D + > D_ and s > r > 1. TTie difference 

U(-±D_,D + )-L(-,\,D_,D + ) 
r 2 s 2 

is monotonically decreasing for x > sexp(1.01-D + -0.01.D_). 

The difference in the latter lemma multiplied by appears if we try to bound 
above ^.{x/r) — fi bi (x/s). Notice that the latter difference is not monotonically 
decreasing from any x onwards, although it can be bounded above by a function 
that is monotonically decreasing for all sufficiently large x. 
Our proof of Lemma ^ uses the following lemma. 

Lemma 5 Let y and 5 be non-negative real numbers. Then the inequality 

+ 1 + KV + 6-2)(y + l) 2 < ^y(y + 3)(y + 5)\ (14) 
holds if either 5 < 2 or y > 0.0099945. 
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Proof. On replacing the inequality sign in fli"^) with the equality sign and squar- 
ing both sides, we obtain an equation of an algebraic curve. Using continuity and 
e.g. Maple's function f solve (for numerically determining roots of polynomial 
equations), the result can then be deduced. □ 

Remark. For y = 0.0099944 and 5 ~ 5.4 inequality ([H]) is not satisfied. Indeed if 
we square both sides of the inequality and take the difference then, considered as a 
polynomial in y, the discriminant has 275 5 - 1985 4 + 4105 3 - 9365 2 + 12995 - 730 
as a factor, which has 5.44694735 • • • as its largest real root. Considered as a 
polynomial in d, we find 

27y 8 - 72y 7 - 2380y 6 - 12792y 5 - 33822y 4 - 48888y 3 - 32076?/ 2 - 2376y + 27 

as a factor of the discriminant, which has 0.00999445028 • • • as its one but second 
largest real root. 

We can now prove Lemma [|. 

Proof of Lemma |]. Differentiating U(x/r, |, LL, D + ) — L(x/s, |, D_, D + ) yields, 
after some tedious calculations, that the derivative is non-positive provided that 
(jg) is satisfied withy = (log(ar/s)-D + )/(D+-D_) and 5 = log(s/r) /(£>+-£>_)■ 
The result then follows on invoking Lemma [5|. □ 



5 Numerical evaluation of certain constants 

For our proof of Theorem |l| we need to evaluate the constants C bl , C b3 , B bl and 
B b3 with enough numerical precision. The purpose of this section is to achieve 
this. (Recall that B b% = \\va x ^(^ n < x h bi {n) / n - (logx)/2).) 

We first consider the evaluation of C b . 3 and C bl (defined by (|TT|) ). We have, 
for Re(s) > 1, 



L fe3 ( S )=(i-3-)- 1 n (l-p-r 1 n 

p=l(mod 3) p=2(mod 3) 



2s\-l 



and 

L bs ( S ) 2 = asMs^^a-s-r 1 n a-^r 1 - (is) 

p=2(mod 3) 

From this, (p]), lim s ^i + o(s — 1)C( S ) = 1 an d the fact that r(|) = y/n, we obtain 

d^j+MA n (1 

p=2(mod 3) 

where for any fundamental discriminant D, xd denotes Kronecker's extension 
(D/n) of the Legendre symbol || Chapter 5]. If x is a real primitive character 
modulo k and x( — 1) = — 1> then 

k 
n=l 



by Dirichlet's celebrated class number formula (cf. equation (17) of |6|, Chapter 
6]). We infer that L(l,%_ 3 ) = w/yTf. Using that Ct s must be positive and 
((2) = 7r 2 /6, we then infer that 



^.=7^71 n (i- P - 2 r" 2 =^i n (i-^ i/2 

V p=2(mod 3) p=l(mod 3) 



Likewise, using that L(l,x_ 4 ) = ir/4, we find formula (|J) for C^. 
Note that, for 9fc(s) > 1/2, 

n n (i-p- 4 -)-'. (i6) 

p=3(mod 4) p=3(mod 4) 

By recursion we then find from @ and (pl|) the following formula, 

l/2 n+1 



which was already known to Ramanujan [|], pp. 60-66] and by Shanks pl| , p. 
78]. Using this expression one computes that C bl = 0.76422365358922066299 ■ • • 
Similarly one can show that 

1 1 °° / _ 2 „ C(2 n ) \ 1/2 ™ +1 

Cb3 = Tis^n (a - 3 ~ 2 \(2", X _ 3 )J ' 

and use it to compute that C& 3 = 0.63890940544534388225 • • • which is in agree- 
ment with the first seven (out of eight) decimals computed for Cb 3 by Shanks and 



Schmid 22 



On noting that, for Re(s) > 1, 



A(n) - 1 C(s) 

L— — = ~cw _C( ) ' 



n=l 



and using that C(s) = l/(s — l)+7 + 0(s — 1), where 7 denotes Euler's constant, 
is the Taylor series for ((s) around s = 1, one infers that 

y*A(n) = £ 1 _2 7 + o(l) = loga;- 7 + o(l). (17) 

Taking the logarithmic derivative of fll~5|) one obtains that 

M?) _ C'(s) L'(s, X - 3 ) log3 v logp 



l 63 ( s ) COO ^,x-3) 3--i p ^ d3)P --r 

from which one easily infers that 

o ^ A fe3 (n) _>^ A(n) L 7 (l,X-3) , log 3 logp , „ m 

^ n -2^ n L(l, X - 3 ) 2 + ^ p 2_ 1 +°l i ^ 

n<x n<x p=2(mod 3) ^ 
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which yields, on invoking (|TT|), 



Similarly we deduce that 

2B bl = -j--j- - + log2 + 2 > — — - 

Note that the argument above yielded that 

L' bi (s) 



B b = - lim -^f4 + -r- 



This can be alternatively deduced from Serre's ]20| proof of @, cf. [11 



As to the numerical evaluation of, for example, the latter prime sum we note 
that 

„ \- logP ^d_ } TT ( 1 \ 

^ P 2 -l ds ° g J-l U-p-2-i a =1 ' 

p=3(mod 4) ^ p=3(mod 4) v r / 

Then, applying ( |i~6"|) m times, we obtain 

£'(2 m ,X- 4 ) C'(2 m ) log 2 

p=3(mod 4) p=3(mod 4) 

Similarly we have 



\ ^ logP \- logp lA[ 

3(mod 41 ^ o=3(mod 4) ^ n=l k 



£(2 m ,X- 4 ) C(2 m ) 



)2 m 



\ ^ logP \- logP I 1 V ( 

p=2(mod 3) p=2(mod 3) n=l v 



L / (2 m , X - 3 ) C(2 m ) log 3 
L(2 m ,x_ 4 ) C(2 m ) 3 2m -: 



Using these expressions, one computes that B bl = 0.163897318634581595856- ■ ■ 
and similarly B ba = 0.1535522449949958272447 • • •. 

Now we can invoke |10|, Theorem 4] to compute the constants C bl (2) and 
C 6s (2). They are given by C f {2) = (1 + B f )/2 for / e {b u b 3 }. We thus find 

C bl {2) = 0.581948659317290797928 • • •, C h {2) = 0.576776122497497913622 • • • 



In [^2] the authors write (in my notation) l B 3 (x) remains so closely propor- 
tional to B\{x) that it is not clear from this data whether C ba (2) > C bl (2) or 
C bl {2) < C bz {2). It would be unlikely that they are exactly equal.' We thus have 
resolved this matter. 

The numerical data from Table 1 in conjunction with the values of C bl , C bl (2) 
and (H) suggest that C bl (3) > and C bl (4) < 0. Similarly it seems plausible that 
C 6s (3) > and C 63 (4) < 0. 



9 



6 Intermezzo: On a claim of Ramanujan 



In the previous section we have seen that B b3 < log a/3. This knowledge suffices to 
disprove a claim that was made in a celebrated, hitherto unpublished, manuscript 
of Ramanujan || on the partition and tau- functions. 

Let r denote Ramanujan's tau-function. Put T n = if 3|r(n) and T n = 1 
otherwise. In Ramanujan's manuscript we read p. 64]: 'We can show by 
transcendental methods that 

y>*=? f n -¥=+o r n 



3 J 1 y/\ogx \(logn) 



where r is any positive number and 

2 l/2 1 _ 7 -2 ! _ 13 -2 X _ lg -2 



3 l/4- ! _ 7-.V ! _ 13 -3- l _ 19 -3 {(! _ 2 -2)(! _ 5-2)(! _ n -2) . . .}l/2 



2, 5, 11, ■ ■ • being primes of the form 3k — 1 and 7, 13, 19, ■ ■ • being primes of the 
form 3k + 1'. This implies that for almost all n, r(n) is divisible by 3. 

Using that r(n) = no\ (n) (mod 3), where o"i(n) denotes the sum of the positive 
divisors of n, it is easy to see that T n is multiplicative and that 



n=l p=2(mod 3) p=l(mod 3) 



From ( |i~9"l) it is not difficult to verify Ramanujan's claim regarding the value of 
C. By logarithmic differentiation we obtain from ( |T9"D that 



y> A r (w) = 21ogp 

n=l p=2(mod 3) p=l(mod 3) 



logp 31ogp 



J9 S + 1 p 



,3s 



On comparing this series with that for Yl^Li ^-b 3 { n ) n s , it is easily seen, on using 
that B b , A < log \/3, that 

b t =b, a - y: + 1) iogp og ^ < ^ < 

indeed we have 5j = —0.53 • • ■. This shows that 

C n / 0.23- ■■ 



VT fc = l + u -^°"" | Q 

3 v/Iogn V logn (log n) 1-1 



where 0.23- ■ • = (1 + Bj-)/2 ^ 0.5 (here we invoked Theorem 4 of |H|) and 
e > 0. Thus the above claim of Ramanujan is false for every r > 3/2 and true 
for r < 3/2. 

Note that if it would be true that Bt = 0, an amazing identity for Euler's 
constant would result. The manuscript H contains several further assertions of 
the type (|i~8l) (with 3 replaced by various other primes), all of which are disproved 
for r > 3/2 in ITJ. 
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7 On the behaviour of 



\( n ) logx 



n<x n 2 

Put Hi(x) = J2 n<x A-bi (n)/n — log \/x, for % — 1 and i = 3. A good understanding 
of the behavior of Hi is needed in order to apply our key lemma, Lemma |[ Let 
us define for i = 1 and i = 3, C+ib/) = sup x>1 iJj(x) and C_(6j) = inf^i H^x). 





Figure 1 Figure 2 

Plot of Hi(x) for 1 < x < 35 Plot of H 3 (x) for 1 < x < 35 



As in [15]| we define V(x;d, a) = n ^ x A(n)/n. It can be shown that as 

n = a(mod d) 



x tends to infinity V(x; d, a) — \ogx/(p(d) tends to a limit C(d, a). Ramare |15| 
has established the following result. 

Theorem 3 [0. For x > 68 we nave \V(x; 3, 1) - ~ logx - C(3, 1)| < 0.1205 
and |V(x;4,l) - | logx - C(4, 1)| < 0.0961. 



We recall from 10 that 



E 



p=a(mod d) 



^<^forx> 289, 



(20) 



and that, for every fixed v > 1 and every x > 

[log X I log D 



i^(l-^)< 

f — 1 X 



E 



r=l 



logf logf 

f r _ f — 1 



(21) 



Theorem 4 VFe /iai>e 

a) C_(&i) = -logV2 and C+(&i) < 0.2663 

b) C_(6 3 ) = -logV3 and C + {b 3 ) < 0.276 

Proof. After some computation for the interval [1, 68] we infer, from Theorem || 
that C+(6i) < B h + 0.0961 and similarly C+(b 3 ) < B fe3 + 0.1205. For the de- 
termination of C~(bi) we use ( p0|) and (|2l| ) in addition to Ramare's inequalities, 
which yields, for x > 289, that Hi(x) > |logx + B bl — 1.3/ y/x — (log4)/x and 
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H 3 {x) >>gx + B b3 - 1.3/ y/x - (log27)/(2x). 



□ 



Let RH(cf) be the hypothesis that for every character x mod d every non-trivial 
zero of L(s, x) is on the critical line. 

Theorem 5 We have 

a) C+(bi) = #i(461) = 0.1701069880305239 ■ ■ •, under RH(A). 
b) C+(b 3 ) = # 3 (3739) = 0.1554480047272349 • ■ under RH(3). 

Proof, (cf. pTOj , Theorem 6]). We recall from [[[(]] that for d < 432 and (a, d) — 1, 
there exists a constant c^q such that for x > 224 we have, on RH(g(), that 



£ 

n = a(mod d) 



A(n) \ogx 
n <p(d) 



< 
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32'Kyfx 



{3 log 2 x + 8 logx + 16}. 



Under RH(4) it follows from this that C + {pi) = max„ a < 6 15xl0 8 Hx(vi), where 2 = 
V\ < Vi < ■ ■ ■ are the consecutive prime powers that can be written as a sum of two 
squares. Similarly under RH(3) we deduce that C + (b 3 ) = max^Lossxio 10 H 3 (wi), 
where 3 = w\ < W2 < ■ ■ ■ are the consecutive prime powers that can be repre- 
sented by the form X 2 + 3Y 2 . On computing these maxima (for details see Section 
0), the proof is then completed. □ 



The reason that, even on GRH, it requires a lot of computation to determine 
C+(bi) and C + (b 3 ) is that these values are so close to B^, respectively Bf, 3 . A 
similar phenomenon occurs in [O for some of the functions considered there (cf. 
Theorem 6). 

Using Theorem ^ and Lemma |2] together with good enough approximations 
for C bl and Cb 3 , one infers that fib L (x) > /i& 3 (x) for x > 27500. After some com- 
putation we then deduce that Conjecture 3 holds true. 

Unfortunately establishing Conjecture 2 requires quite a bit more work. In 
particular we need values for D_ and D + in Lemma that are more closely to- 
gether than those coming from Theorem |4j. Without improvement of Theorem |3|, 
the upper bounds in Theorem ^ cannot be improved. The lower bounds, however, 
are amenable to improvement. 

Let Af(x) denote the quantity that is sandwiched between D-fif(x) and 
D + fif(x) in (|9|). Using the lower bound for H 3 (x) appearing in the proof of 
Theorem £|, we deduce that H 3 (x) > for x > 25. We infer that 



A ba (x) > -\ogV3{ 



fx b3 (x) 



A**(^)}. 



On applying Lemma || with D_ = —log a/3 and D + = 0.276, we deduce that 
A b3 (x) I n H {x) > -0.09586- •• for x > 10 9 . Taking D_ = -0.09586 ■• • as new 
value and repeating the procedure, we obtain D_ = —0.06890 ■ • -. Iterating twice 
more, we see that for x > 10 9 we can take _D_ = —0.0672 in Lemma |^. 

For any x satisfying the conditions of Lemma |2], we can proceed as above. 
If the first iteration yields an improved value of our initial D_ (which we take 
to be — log a/3), then it is not difficult to see that every further iteration yields 
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a value of _D_ not less than the previous one (this is so since, for given r > 1, 
L(x/r, ~, D_, D + )/U(x, |, D_, D + ) is increasing, considered as a function -D_). 
On the other hand the value cannot be improved beyond zero and hence the 
iteration process must converge. If the first iteration does not yield an improved 
value for D- (which is initially taken as — log V$), we put Wi(x) = — log a/3 for 
every i > 0, otherwise we put Wq(x) = — log v3 and define 



,, . .. . H%\M& 0-276) \log3 
!< 1 r(,-.i. /r,(.r). ().27(i) 



Empirically it seems that after n iterations we can expect to have approached 
the limit value lim^oo w%(x) with 0(n) decimal precision. 

For bi we proceed similarly. After some computation using the lower bound 
for Hi(x) given in Theorem |], we find that Hi(x) > 0.065 for x > 97. Hence 

H x {x) > (- log v^- 0.065) {fi bl {x) - fi bl {^)} + 0.0Q5fi bl {x). 

If the first iteration does not yield an improved value for D_ (which is initially 
taken as — log a/2), we put Vi(x) = — log a/2 for every i > 0, otherwise we put 
vo(x) = — log v2 and define 

, . fL(-§f.\,vAx\ 0.2663) \ r 

d ^ = tt( ~ n J ~ 1 (1 ° g ^ + °- 065 ) + °- 065 - 
\£/(x, jj,Vi{x), 0.2663) J 

To sum up, we have established: 

Lemma 6 Suppose that x > xo > 2 and i > 0. T/ien @ /io/ds irue wn£/i f = b±, 
D- = Vi(xo), D + = 0.276. It also holds true with f — 63, £)_ = ?Z>j(xo) and 
D + = 0.2663. 

This lemma, although amenable to further improvement, is sufficiently sharp for 
our purposes. 

8 The proof of Theorem @ 

Before proving Theorem 0, we will need two more lemmas. From prime number 
theory we recall that i/j(x; d, a) = J2 n <x, „E fl (mod d) A W- 

Lemma 7 We have 

a) Vfci 0*0 > 0.4924x /or x > 37. 

b) ip b3 (x) < 0.5176a; for x > 3793. 

Proof. Let d < 13 and (a,d) = 1. Then |^(x;<i, a) — x/(p(d)\ < \/x for 
224 < a; < 10 10 by gS], Theorem 1] and |^(x;d,a) - < 0.004560x/<^(d) 

for x > 10 10 by [16, Theorem 5.2.1]. From these inequalities the lemma follows 



after some computation. □ 

For y > 3 we define Sb 3 (y) by 0.5176?/, except for the intervals [3,49), [49, 181), 
[181, 487), [487, 1369), [1699, 1933), [2287.2437), respectively [3733,3793), where 
we define S bs (y) to be (respectively) 0.653954?/, 0.605778?/, 0.557372?/, 0.534528?/, 
0.526579?/, 0.521825?/ and 0.51996?/. 
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Lemma 8 For y > 2 we have ipb 3 {y) < Sb 3 {y)- 

Proof. The points where tpb 3 and S& 3 change value occur only at prime powers 
representable by X 2 + 3Y 2 , which we denoted by 3 = W\ < w 2 < ■ ■ ■■ We now 
check that ipb 3 { w t) — Sb 3 (wi) for every Wi < 3793. For Wi > 3793 the result 
follows by Lemma |7|. □ 




20 40 60 80 100 



Figure 3 
Plot of ip bl (x) versus 0.4924x 
for < x < 100 




Figure 4 
Plot of ipb 3 {x) versus Sb 3 (x) 
for < x < 100 



At last we are in the position to prove Theorem |2]. 

Proof of Theorem |2|. As we have shown in Section 2, it suffices to establish 
Conjecture 2. To this end we have to prove that, for x > 8, 



Let us denote the 6 intervals in the definition of Sb 3 (y) by [r^, Sj) for i = 1, • • • , 6 
and put CKj = S^ir^/r^ — 0.5176 (note that a, > 0). From Lemma § we infer that 



K( x ) < z2 b z(n)S b3 (-) = 0.5176/^ 3 (-) + VaiK(-) - Vb 3 (— )} 



i=l 



Si 



Put £0 = 1.5 x 10 11 . Using a computer (see Section |^) Conjecture 2 can be es- 
tablished for x < xq. Hence now assume that x > xq. For notational convenience 
we shorten U(x/r, ^,Ws(xo/r), 0.276) to U 3 {x/r), L(x/r, |, wg(xo/r), 0.276) to 
L 3 (x/r) and L(x/r, ^,Vs(xo/r), 0.2663) to L\(x/r), where r is some fixed number. 
On applying Lemma || we deduce that 

X b3 (x) 



2C> 



b 3 



< 0.5176[/ 3 (^) + V Oi{l/s(-) - L 3 (-)} 



i=l 



By Lemma § each of the six terms in the above sum is non-increasing for x > xq 
and thus the sum is bounded above by its value in xq, which on its turn is less 



14 



than 0.0224£/3(a;o/3). One easily checks that U${x/$) > U^{xq/3) for x > xq 
(on noting that U(y,r, D-, D + ), considered as a function of r, is increasing for 
y > exp(L> + + (D + - D_)/t)). We thus obtain that \ bz {x) < \mC hz U z {x/3). 
Using Lemma |6| and the lower bound for ip^ given in Lemma [j], we infer that 
> EnKx/sjbM^ix/n) > 0.4924/i 6l (x/37) > 0.9848C 6l Li (x/37). A 
computation shows that 0.9848(7^1 (xq/37) > 1.08Cb 3 U3(xo/3). By Lemma |3| 
we then have 0.9848C bl L bl (x/37) > 1.08C h U 3 (x/3). for every x > x . We thus 
obtain that for every x > Xq, 

X h (x) > 0.9848C 6l L x (^) > 1.08C ba U 3 {^) > 
completing the proof. □ 



9 Computations of results used in Theorems g 
and [5] 

In the proof of Theorem we have used that Conjecture 2 is true for x < xq with 
xq = 1.5 x 10 11 . We have established that result as follows. 

Checking Conjecture 2 requires the computation and comparison of the sums 

= ^2bi(n) log n, i = l,3, 

n<x 

and, consequently, the computation of the characteristic functions b\{n) and 63(71) 
for all positive integers n < xq. Because of the size of Xq, the range of a;- values 
for which Conjecture 2 had to be checked, was split up in sub-intervals of length 
10 6 , large enough for efficiency, and small enough to avoid so-called cache misses 
during the computations. 

We first describe the case A&^n). For a given interval, say, [A, B], an integer 
array b(j),j = 1, 2, . . . , B — A + 1 of length B — A + 1 is initialized to 0. Here, 
b(j) corresponds to b±(j + A — 1). Next, all the possible sums of squares x 2 + y 2 



y e 



^A/2,VB 



of integers < x < y, with A < x 2 + y 2 < B, hence x G 0, yB~/2 

are computed as follows. First, the sequence of all the squares 

y 2 G [A/2, B] is precomputed and stored. Next, for each x — 0, 1, . . 
the sums x 2 + y 2 are computed 

for all y 2 e [max(A - x 2 , A/2), B - x 2 ] C [A/2, B]. 



VBj2 



For all the sums x 2 + y 2 =: n obtained in this way, b(n — A + 1) is set to 1. 

The case Xb 3 {n) is treated similarly: the same initialization of array 6 is carried 
out. Next, all the possible sums x 2 + 3y 2 of integers x, y, with A < x 2 + 3y 2 < B, 

hence x G 0, \f~B , y E 0, \[B~j3 are computed as follows. First, the sequence 

of all the triples of squares 3y 2 G [0, B] is precomputed and stored. Next, for 



each x = 0, 1, 



the sums x 2 + 3y 2 are computed 



for all 3y 2 G [max(0, A - x 2 ), B - x 2 ] C [0, B] 
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z/10 
1 

1.5 
2 
3 
4 
5 
6 
7 
8 
9 

10 



A&i(g) 

378458908590.818 
572353849423.260 
767521856517.400 
1160486988190.213 
1555965223692.576 
1953301629004.525 
2352112868630.901 
2752146230205.959 
3153223047545.408 
3555209733889.339 
3958003171956.632 



K( x ) 

316358774044.179 
478438468735.511 
641582406621.494 
970068358550.987 
1300655152892.098 
1632795521743.015 
1966168966371.294 
2300563843364.554 
2635831188875.970 
2971859287714.156 
3308561817015.470 



B^x) 
15570512744 
23160971166 
30700929088 
45678037444 
60558145064 
75367348255 
90120785046 
104828319151 
119496904413 
134131682979 
148736628858 



B 3 (x) 
13015595425 
19360573686 
25663340448 
38182949191 
50621477125 
63000746043 
75333407591 
87627692348 
99889427349 
112122909167 
124331455166 



Table 2: X bl (x),Bi(x) versus Xb 3 {x), B 3 (x) 



and for all the sums x 2 + 3y 2 =: n obtained in this way, b(n — A + 1) is set to 1. 
This corresponds to b 3 (n). 

We have implemented these algorithms for bi(n) and 63(71) in Fortran and 
used them to compute X^x) for i = 1,3, and to verify Conjecture 2 for x = 
8, 9, ... , 1.5 x 10 11 on one 250 MHZ processor of CWI's SGI Origin 2000 computing 
system. Computing time was 7.6 CPU hours. We also used our program to check 
the values of B\(x) = ^2 n < x bi(n), gi ven f° r x — 10% i = 1, . . . , 12, by Shiu in 
Table 1 of P31 (where B\(x) is called W(x)). Computing time to extend our 
results from 1.5 x 10 11 to 10 12 was 77 CPU hours. We found agreement with Shiu 
for i = 1, . . . , 10, but differences for i = 11 and % = 12: ^(lO 11 ) = 15 570 512 744 
and ^(lO 12 ) = 148 736 628 858, whereas Shiu gave W(10 u ) = 15 570 523 346 
and H^(10 12 ) = 148 736 629 005. Shiu used a different, more efficient method than 
ours, but he has confirmed our value of _Bi(10 n ) after checking and correcting 



his program |24 



We have spot-checked our program for computing bi(n) and 63 (n) on various 
intervals of length 10 6 with the help of Lemma |]. This requires the decomposition 
in primes of each n for which we wish to compute b\ (n), which is extremely 
expensive, compared with composing all integers in a given long interval [A, B] 
as a sum of integer squares. However, we found agreement for all the checks 
we did, in particular for those in the neighboorhood of x = 10 12 . In Table 2, 
we list, for i = 1,3, the values we found of X^x) and B^x) for x = j x 10 11 , 
j = l,1.5,2,--.,10. 

In the proof of Theorems 5a and 5b, we have used that 

max HAvi) = #i(461) = 0.170106..., , respectively, (22) 

t) 4 <6.15xl0 8 

max H 3 (wi) = # 3 (3739) = 0.155448..., (23) 

Wi<1.083xl0 10 

We have established these results as follows. 

Let x = 6.15 x 10 s . We first generated the primes < \fx with the sieve of 
Eratosthenes, and stored the following pairs (n, A^n)): 

(2 fc ,log2),fc = l,2,...,Llog 2 a;J, 
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(p 2k ,2\ogp),k = l,2,... 



2 1o Sp x 



for the primes p = 3 mod 4 < \/x, 



(p k , logp), k = 1, 2, ... , [logp xj , for the primes p = 1 mod 4 < \/ir, 

into an array, sorted increasingly according to the first element of the pairs. 
The set of numbers n in these pairs in fact contains as a subset all the prime 
powers f i, t>2, . . . < \fx which can be written as a sum of two squares. For these 
(n, Aft 1 (n))-pairs, we computed Hi(n) and we verified that 

max Hxfa) = #i(461) = 0.170106... 

The remaining interval [[V6.15 x 10 8 J + l,x] was split up in pieces of length 10 , 
and for each of these intervals, [A, B], say, the primes p = 1 mod 4 were generated 
with the sieve of Eratosthenes, together with logp. These pairs (p, logp) were 
mixed with the (n, A& 1 (n))-pairs generated above for which n G [-A, S] and then 
it was verified that max^^^^j Hi(vi) < Hi(A61). This proved (^). Computing 
time was 81 CPU seconds. Relation ( p3[ ) was proved in a similar way at the 
expense of 1340 CPU seconds. 



10 An alternative approach 

In the previous sections we have made essential use of asymptotic information 
regarding the distribution of primes. Some of the results we used depend even- 
tually on RH(3) and RH(4) to be true up to some finite height. It might come 
as a surprise then that it is possible to show that B\(x) > B^{x) for x > 10 9111 , 
without invoking any result from computational prime number theory (one only 
needs the ability to compute some successive primes...). 

Our method of establishing this is inspired on Selberg's [|1^, pp. 183-185] 
method of obtaining an asymptotic evaluation for N(x;4, 1), where N(x;d,a) 
denotes the number of integers n < x that have no prime factor p with p ^ 
a(mod d). Unfortunately Selberg's method does not seem to generalise well; for 
example we have no idea how to generalise it so as to show that N(x; 4, 3) > 
N(x; 4, 1) for x > Xo, with xq some effectively computable constant. Vide (TT 
for generalisations of Selberg's method. 

Lemma 9 a) For x > 2 we have 

IB^x) -C bl —^=\ < 9.62- " 



^/\ogX hgx' 



b) For x > 2 we have 



\B 3 (x) - C 63 X I < 8.53; 3 



Vlogx log X ' 



Corollary 1 For x > 10 9111 we have Bi(x) > B 3 (x). 

In the proof of Lemma we will make use of the following result. 
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Lemma 10 Let C2 = 2e 7 and C3 = V3e J . For z > 1 we put 

1 z 1 2 

f( z ) = z z2 o lo g( c 2^) and 9{z) = z 22 -zlog(c 3 z). 



n 2 ^— ' n 3 

n<2, 2fn i<£, 3fn 



sup \f(z) \ = -/(3 - ) = ^{log6 + 7} - 3 = 0.55346270119438- •• 

2>1 2 

and 

o 

sup |#(;z)| = -0(4-) = -log(4v / 3e 7 ) - 6 = 0.70084312094794 ■ • • 

z>i 3 

Proof. We only prove the statement concerning f(z), the statement regarding 
g(z) can be proved in a similar way. 

Since z log(c2z) is monotonically increasing, we obtain, cf. [|X0j, Lemma 4], that 
sup^l/M = sup{|/(l)|,|/(3-)|, |/(3)|, |/(5")|, |/(5)|, •••}. Using the Euler- 
MacLaurin summation formula, cf. [^, p. 6], one finds that for integers n > 1, 

1 11 0i(n) , > 

5:»= iog " +7+ ^-i2^ + 6^- (24) 

m<n 

where 0i(n) G [0,1]. Clearly 

y - = y--- y (25) 

m<n, 2\ra m<n m<n/2 

Let n > 3 be an odd integer. Notice that f(n~) = f(n) — 1. Using (|24"D and (|25"D 
it is not difficult to deduce that 

1 x 7i , / n \ 7i 1 
-- < /(n) < - log( + + 



2- JW -2 n-L 6(n-l) 2 60n 3 
and 

1 ,, x 1 1 1 2 

— < -f( n ~) < - H 1 1 . 

2" M ' ~ 2 2(n - 1) 12n 2 15(n - l) 4 

Using that the latter two right hand sides are monotonically decreasing in n and 
noting that sup^ \f(z)\ > |/(3 - )|, we see that sup 2 > x |/(*)| = sup^^n \f(z)\ = 
l/(3")|. □ 



Let {z} — z — [z] denote the fractional part of z. Using (Ell) it can be shown 
that the functions / and g are almost periodic in the sense that they converge 
uniformly to the periodic functions ~ — l 5 ^-}, respectively 1 — {^-} — {^-}, 
(cf. Figure 5). 
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Figure 5 Figure 6 

Plot of g(x) for 1 < x < 10 Plot of g(x) minus its limit function 

for 1 < x < 20 



11 Proof of Lemma g 

Let P2,P-s denote the set of primes p that satisfy p = 2(mod 3), respectively 
p = 3 (mod 4). Let (P?), (-P3) denote the set of natural numbers that have no 
prime divisor p with p ^ 2(mod 3), respectively p ^ 3(mod 4). Let ip 3 (x),ip4:(x) 
denote the number of integers 1 < n < x that have no prime divisor p with 
p ^ 2 (mod 3), respectively p ^ 3 (mod 4). 



Proof of part a. Put c 2 = 2e 7 . We consider the expression 

00 

EE E E^)io g ^ = 

j=Q m£(P3) l<n<x/(23m 2 ) d\n 

n = l(mod 4) d £(P 3 ) 

OO 

EE E M**£r d E i- (»> 

J=0 me(P3) de(P 3 ) d | n> i<„< I /(2i m 2) 

d<a;/(2im 2 ) n = l(mod 4) 

By approximating both sides of this equation in terms of the function B\, we will 
arrive at an approximate functional equation, (P4"|), for B\ which on solving will 
yield an explicit lower bound for B\. 

On recalling that ^ d \ n n(d) logrf = — A(n), one sees that, when n = l(mod 4) 
and d < x/(2- J m 2 ), we have 

flog ^2% ifne(Pi); 

w(fi) log — — — = < logp if n is divisible by exactly one p from P 3 ; 
'2P tti d I 

d|« " 10 if n has > 2 distinct prime factors from P3. 

On noting that 



£i(*) = E E ^ 
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3=0 me(P 3 ) 
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we see that the left hand side of equals 



oo 

i^logc^+^logpS^-E E M2 J m> 4 (^), (27) 

peP 3 P j=0 me(P 3 ) 

?'>1 

which we write as 

Bi (x) log c 2 x + / x (x) - I 2 (x) . (28) 

We write z = x/2^ and consider the expression formed by the three inner sums 
in fl26|) , that is 

E E /•m**^ E i- (») 

mS(P3) d e( p 3) d|n, l<n<z/m 2 

d<z/m? n = l(mod 4) 

Given an integer /c, let a"o(fc) denote the product of the distinct primes that occur 
to an odd power in the prime factorisation of k. We put <Jo{k) = 1 if if there is 
no prime that occurs to an odd power in k. Note that 

E M = ri°o(k)), (30) 

m 2 d=k 

where the sum is over all integers m and d such that m 2 d = k. On writing 
m 2 d = k in (p9|), and invoking (0), we deduce that the triple sum in ( p9|) equals 

EM^))io g ^ ^ l 

fce(p 3 ) ' kx<z/k 

k<z fc 1 =fc(mod 4) 

The right hand side of ( PEP is thus seen to equal 

oo 

E E M-oW)iog|f e i. 

j=0 fe6(P 3 ) fe 1 <i/(2J'fc) 
k<x/23 ^Etlmod 4) 



which simplifies to 



E »{<d'))\og C -f 

de(pp k 1 <x/d 

d<x fc^i'lmod 4) 

where d' denotes the largest odd divisor of d and Pg = P 3 U{2} and (Pg) is defined 
as (P3) but where now no prime divisor p with p = l(mod 4) is allowed. The 
right hand side of (p6|) is thus seen to equal 

E 



d'-i 



x 2+ -1 — 
1 — 

Ad 4 



^(«T (ti))log— = - 2^ ^ log— + / 3 (a;) 



d<x d<x 

(31) 

where 

|/3(x)| < \ E tog^- (32) 



20 



For 1 < d < x we have, recalling the definition of / (made in Lemma |i0|) , 

1 , c 2 x 1 i- f x ^d 

2 lo ^= E n'^x 

n<x/d, 2-fn 

Combining the latter equation with the sum in the right hand side of ( |5T| ) yields 
x sr^ fj,(a (d')) c 2 x x ^ n(a (d')) , . 



4 ^ d d 2 ^ 



where 



Note that 



< \ J2 1/(3)1- ( 33 ) 



2 ^ u v d 

d<a; 



dn<x, dG(Pg) fc<X dn = k, dG(Pg) 

2fn 2-fn 

Denote the latter inner sum by h(k). We claim that = b\(k). First let us 
consider the case where k is odd. Then 



p-i 

! 2 



dn=k, de(P3) dn=fc p|d 

where u(d) denotes the number of distinct primes dividing d. We see that for 
odd k, h is the Dirichlet convolution of two multiplicative functions and is thus 
itself a multiplicative function. For arbitrary k we note that h(k) = h(k'), where 
k' is the largest odd divisor of k. Thus h is a multiplicative function. An easy 
computation shows that for every prime power q we have h(q) = b\(q). Since 
both h and b are multiplicative, this completes the proof of the claim. We thus 
infer that 

x_ ^ MgoW) = xy* Mm) = x r dgift) = : Bi(x) x r gift) 
2 ^ 2 ^ m 2 L t 2 2 /, t 2 

2-fn 

Thus the right hand side of (^) equals 

~2 + 2 / + 3{X) + h{x) ' 



(34) 



Equating it with the expression ( Bq ) for the left hand side of (Eq), we get 
B^x) hgc 2 x ^-dt = -I x {x) + I 2 (x) + h{x) + h{x) + 



Next we will consider effective estimates for Ij(x) for 1 < j < 4. Using the 
trivial estimate B\(x) < x, we obtain that 

< h(x) < x V < .23a:. 
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On noting that 



E E 



log^m 2 



2 E 



log(2m 2 



< 2.7 



j=0 me(P 3 ) me(P 3 ) 

and ipi(x) < x, we deduce that < h^x) < 2.7 'x. Using ( |32"1) we deduce that 

i t / \ i 3 v — a , cx 3 v — a f dt 3 _ __ 
|Js(x)|<- lo §7^E y < iC x<2.68x. 

<fc(P'), rf<x i<d<x Jd 



For J^x) we have, by (|33|) and Lemma [T0|, | _Z" 4 (a;) | < 0.277x. 

Put = f* B^dt/t 2 . An easy calculation (divide by x 2 log 3//2 x and 

integrate) now shows that if 



x 



— a_x < x \ogxA!(x) A(x) < a + x, for x > x 



(35) 



then there exists a constant Co such that 

c a/ log x — 2a + < A(x) < c ^/\ogx + 2a_, for x > x . 
On inserting the latter estimate in ([35|) and invoking ([!]), it then follows that 



Pi(x)-C 6l 



Vlogx 



< (ot_ + a_ 



.r 



log x 



for x > Xq. 



(36) 



From our estimates for Ij(x) with j — 1, • • • , 4, we see that we can take a_ = 3.96, 
ct + = 5.66 and x = 2. □ 

Proof of part b. Making the obvious modifications in the proof of part a, we 
deduce that 



x f x B 3 (t) 
B 3 (x) logc 3 x — — I — - — dt 



■Ji(x) + J 2 (x) + J 3 (x) + Ji(x) + 



(37) 



where Ji(x) = E pe p 2 , r>i lo gP 5 s(^f), ^(z) = E^o E me( P 2 ) MS^ 2 )^^?), 



E ^ 



>C and \J A (x)\ < ]- ^ \g( X 



de(P^), d<x 



d 



2 ^ [JK d J 
de{P£), d<x 



with P' 2 = P2 U {3} and (Pj) defined as (P2), but where now no prime divisor p 
with p = l(mod 3) is allowed. Reasoning as before we find 

< JAx) < x V < 0.36x and < J 2 (x) < ^-x V log ( 3 ™ ) < 2.7x. 

peP 2 me(P 2 ) 

Furthermore we find that | «/ 3 (a;) | < 2xc 3 /3 < 2.06x and |J,4(x)| < 0.36x. From 
these estimates and fl3?D we infer that we can take a + = 5.12, a_ = 3.41 and 
xq = 2 in the analogue of (|36|) . □ 



22 



Remark. In the proof of part a we have used the trivial estimates ipi{x) < x and 
Bx(x) < x. Using that the integers n counted by ^(x) satisfy n = l(mod 4) and 
3 \ n, we obtain the sharper estimate 



, . . r x + 11., r x + 7. x 7 
— 12 J 1 12 J - 6 6 



(38) 



Similarly, some computation yields that B\(x) < x/2 + 2. In this way the value 
9.62 appearing in Lemma |9| a can still be further decreased, but we have not 
carried this out. Similarly the estimates in the proof of part b can be improved. 
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